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1.  SUMMARY 


The  primary  goal  of  this  project  was  to  compare  three  quantitative  approaches  to  modeling 
EM  induction  for  UXO  discrimination:  a  phenomenological  dipole  model,  a  semianalytic 
theory,  and  a  finite-element  numerical  method.  The  secondary  goal  was  to  quantify  the  value  of 
multiple  spatial  components  and  time  channels.  A  general,  time-dependent  triaxial  dipole  model 
was  developed  that  was  very  successful  in  discrimination  tests.  The  semianalytic  theory  was  not 
completed  and  implementation  of  the  numerical  method  was  impractical  given  the  allocated 
resources.  The  relative  merits  of  additional  temporal  and  spatial  information  were  assessed  by 
analyzing  decimated  data  from  the  Geonics  EM61-3D.  The  responses  of  25  unique  ordnance- 
and  scrap-like  objects  buried  in  45  depth-orientation  states  were  modeled  as  triaxial  time- 
dependent  dipoles  and  were  classified  using  the  properties  of  the  inferred  model  parameters. 
Performance  was  measured  by  the  fraction  of  false  positives  for  ordnance-like  objects  at  91% 
true  positives.  The  best  false-positive  rate  for  the  full  multicomponent,  multichannel  data  was 
just  5%,  which  was  achieved  using  the  power-law  slope  of  the  time  decay  as  a  discriminant: 
ordnance-like  items  remain  in  “intermediate”  time  longer  and  display  a  characteristic  -3/2  slope, 
whereas  scrap-like  items  (even  those  made  of  steel)  transition  to  “late”  time  and  the  recovered 
power-law  indices  vary  widely.  However,  this  excellent  performance  is  largely  due  to  the  fact 
that  the  ordnance-like  objects  are  larger  than  the  scrap-like  objects  in  this  data  set.  Better 
generalization  may  be  obtained  with  discriminants  based  solely  on  shape — ordnance-like  objects 
are  usually  axisymmetric  about  a  long  axis — which  yielded  32%  false  positives,  again  in  the 
power-law  slope.  In  contrast,  the  best  false-positive  rate  for  single-component,  single-channel 
data  was  68%.  Other  instrument  configurations  (single-component,  multichannel  or 
multicomponent,  single  channel)  were  generally  ranked  intermediately  between  these  extremes. 
This  work  demonstrated  in  an  internally  consistent  way  the  high  perfonnance  in  UXO 
discrimination  that  can  be  achieved  with  multicomponent,  multichannel  electromagnetic  sensors, 
as  well  as  the  value  of  relatively  simple  modeling  and  discrimination  procedures. 
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2.  INTRODUCTION 


Electromagnetic  (EM)  induction  sensors  have  been  recognized  as  a  principal  means  of  UXO 
detection  and  discrimination  (e.g.,  Bell  et  a/.,  2001).  Two  issues  regarding  the  use  of  such 
sensors  are  the  kind  of  data  acquired  and  the  algorithms  used  to  interpret  the  data.  Time-domain 
EM  (TDEM),  or  pulsed-induction,  sensors  use  an  abrupt  change  in  transmitter-coil  current  to 
induce  eddy  currents  in  a  target,  whose  decay  is  sensed  in  a  receiver,  usually  another  coil.  Most 
such  instruments  use  a  transmitter  coil  in  the  horizontal  plane.  They  differ  in  the  amount  and 
complexity  of  data,  specifically,  the  number  of  receiver  spatial  components  and  the  number, 
location,  and  width  of  time  gates  that  record  the  eddy-current  decay.  Physics-based  models  are 
used  to  interpret  these  data,  and  three  such  models  were  studied  here.  A  semi-empirical  method 
treats  targets  as  a  group  of  infinitesimal,  orthogonal,  time-dependent  dipoles  and  constructs 
responses  for  arbitrary  position  and  orientation  through  linear  superposition.  Classification  can 
be  improved  by  experimentally  determining  the  directional  responses  of  targets  of  interest.  The 
technique  is  very  fast  but  is  limited  to  distances  relatively  far  away  from  the  source  and  receiver, 
assumes  no  environmental  effects,  and  can  require  prior  data  for  each  object  to  be  classified.  A 
semi-analytic  theory  is  relatively  fast,  models  the  full  field  at  arbitrary  distances,  and  can  include 
ground  conductivity,  but  restricts  target  shapes  to  solid  triaxial  ellipsoids.  Numerical  models 
such  as  the  finite-element  method  offer  the  most  accurate  solutions  for  arbitrary  objects  and 
environments  but  are  slow,  unsuited  to  parameter  estimation,  and  require  that  responses  from  all 
potential  targets  be  cataloged.  The  purpose  of  this  investigation  was  first  to  assess  the  relative 
merits  of  these  three  models  and  second  to  quantify  the  value  of  increasing  the  amount  and  kind 
of  TDEM  data  for  UXO  discrimination,  providing  metrics  for  such  discrimination.  Our 
approach  focuses  first  on  model  development  and  second  on  test-plot  measurements  of  inert 
ordnance,  ordnance  simulants,  and  non-ordnance. 
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3.  MODELING 


Three  different  modeling  techniques  were  proposed  for  this  work.  However,  an  extension  of 
the  simple,  widely  applied  dipole  model  proved  very  effective  for  our  discrimination  study  while 
the  other  two  reached  dead  ends,  at  least  as  far  as  we  were  able  to  pursue  them  here. 

3.1  Mean-Field  Theory 


Maxwell’s  equations  for  the  electromagnetic  fields  in  a  compact,  conductive  object  can  be 
reduced  to  the  single  equation 


Vx 


-Vxew 

.M 


=  a„ae 


.(«) 


(3.1) 


where  a  and  //  are  the  conductivity  and  permeability,  respectively.  This  is  an  eigenvalue 
problem,  where  the  eigenvalues  a„  are  the  mode  constants  and  the  eigenvectors  e"'  are  the 
electric-field  shapes  for  mode  n.  This  electromagnetic  problem  is  therefore  analogous  to  the 
mechanical  problem  of  the  motion  of  a  drum  head:  depending  on  where  the  drum  is  struck, 
different  vibrational  shapes  are  excited.  Similarly,  different  modes  of  electric-field  circulation  in 
the  target  will  be  excited  depending  on  shape,  position,  and  orientation  with  respect  to  the 
transmitter,  and  these  fields  will  be  sensed  differently  depending  on  the  receiver  configuration(s). 
The  time-domain  solution  to  Equation  (3.1)  for  the  electric  field  is 

E(x,t)=  24,eWWexp(-«»0  (3-2) 

n 

where  the  coefficients  A„  give  the  amplitudes  of  each  mode. 

The  essence  of  the  “Mean-Field  Theory”  (MFT),  developed  by  Dr.  Pieter  Weichman  while 
he  was  at  Blackhawk  Geoservices,  is  to  find  a  truncated  Taylor-series  solution  to  Equation  (3.2) 
which  yields  an  approximation  of  the  shape  functions  e(n)  and  time  constants  an.  These  can  be 
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precomputed  based  on  the  target  shape  and  the  model  execution  speed  increased  by  referring  to 
these  parameters  in  a  look-up  table.  An  actual  modeling  run  then  needs  only  to  calculate  the  A„ 
parameters  based  on  the  source-target-receiver  geometry.  A  triaxial  ellipsoid  was  adopted  as  a 
canonical  shape  that  was  analytically  tractable  to  bound  the  shape  polymomials  but  flexible 
enough  to  represent  diverse  objects  such  as  cylinders,  discs,  and  spheres  (Figure  3,1). 

While  elegant  in  principle,  the  MFT  never  reached  practical  application  while  the  author  was 
at  Blackhawk.  Although  the  lowest-order  (“late-time”)  constant  ai  could  be  accurately 
computed  and  compared  to  analytic  models,  computation  time  and  accuracy  limited  how  well  the 
higher-order  time  constants  could  be  determined.  The  author  was  considering  a  one-dimensional 
approximation  that  may  have  provided  a  “patch”  but  would  have  been  inconsistent  with  the 
overall  approach  of  the  model.  Extension  to  penneable  (ferrous)  and  hollow  objects  is  essential 
for  treating  the  UXO  problem  but  was  never  within  sight.  After  the  author’s  departure,  we 
transferred  development  of  the  code  to  collaborator  Alex  Becker  and  his  staff  at  UCB.  This 
work  was  pursued  under  different  funding  and  resulted  in  an  order-of-magnitude  improvement  in 
the  early-time  accuracy  of  the  MFT  but  still  has  not  reached  objects  with  extreme  shapes  or  those 
that  are  steel  and/or  hollow.  While  MFT  development  continues  externally  at  lower  priority,  we 
have  elected  to  pursue  other  alternatives. 
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Figure  3.1  Example  of  fundamental  (n=l)  mode  shapes  for  a  1.25:1:1  nonpermeable  solid  ellipsoid  calculated  using 
4th  order  MFT. 
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3.2  Numerical  (Finite-Element)  Model 


Piecewise  solution  of  Maxwell’s  Equations  leads  to  the  finite-element  method,  which  is 
widely  used  for  physical  analysis  of  complexly  shaped  systems.  We  tested  the  ability  of  the  low- 
frequency  electromagnetic  module  of  the  ANSYS  package  to  solve  UXO  problems.  The  original 
motivation  was  to  test  and  calibrate  the  MFT  by  examining  the  TDEM  response  of  ellipsoids. 
Axisymmetric  equal-volume  objects  with  aspect  ratios  of  1:1  (sphere),  4:1  (prolate  ellipsoid)  and 
1:16  (oblate  ellipsoid)  were  modeled  at  1  m  directly  below  an  EM-61  in  longitudinal  and 
transverse  orientations.  In  the  former,  the  axis  of  symmetry  is  vertical  and  in  the  latter  the 
symmetry  axis  is  horizontal.  Runs  were  performed  for  both  aluminum  and  steel  targets. 
Aluminum  objects  had  late-time  constants  ai  measured  from  the  decay  curves  in  good  agreement 
with  theory.  Steel  objects  are  dominated  by  magnetic  polarizability  and  the  rules  developed  for 
eddy  currents  do  not  apply.  Instead,  the  magnitude  of  the  response  varies  directly  with  the 
amount  of  steel  in-line  with  the  primary  field  (Figure  3.2). 

The  full  capabilities  of  ANSYS  were  tested  by  constructing  a  geometrically  accurate  model  of 
an  81 -mm  mortar  projectile,  consisting  of  a  steel  body,  an  aluminum  fuze  and  tail  assembly,  and 
a  “hollow”  interior  .  The  model  was  evaluated  in  nose-down  and  horizontal  orientations  beneath 
an  EM-61.  The  calculations  reveal  the  detailed  structure  and  evolution  of  eddy  currents  (e.g., 
Figure  3.3)  and  magnetic  induction  in  the  target,  providing  substantial  insight  into  the  physics  of 
EM  sensing  of  UXO. 

The  ANSYS  calculations  provided  the  baseline  for  the  MFT  as  required.  Another  objective 
was  to  be  able  to  generate  realistic  “synthetic”  data,  saving  the  labor  and  data-reduction  costs  of 
experimental  studies — this  is  a  general  goal  of  such  modeling  in  industry.  We  found  two 
obstacles  to  this.  First,  we  discovered  that  working  with  ANSYS  is  a  full-time  job  (perhaps  well- 
known  by  the  consultants  on  whom  we  relied  for  the  early  parts  of  this  work)  and  we  were 
unable  to  devote  the  resources  necessary  to  rapidly  and  robustly  generate  target  “CAD”  models 
and  finite-element  meshes.  Although  there  are  many  tools  available  in  ANSYS,  they  demand  a 
very  large  time  investment  to  leam  and  master.  Second,  the  cost-effectiveness  trade  point 
seemed  out  of  reach.  Model  execution  times  are  very  slow  (hours  for  a  single  position  of  the 
sensor  with  respect  to  the  target)  that  several  days  of  CPU  time  would  likely  be  expended  for 
measurements  that  can  take  as  little  as  a  few  minutes  for  simple  sensors  and  perhaps  an  hour  for 
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a  careful  grid  with  a  complex  sensor.  We  will  reserve  future  finite-element  analysis  for 
understanding  unusual  target  responses  that  cannot  be  adequately  modeled  with  simpler  methods. 


FEA  Results  for  Steel  Sphere,  Proiote  ana  Oblate  Spheroids 


Figure  3.2.  ANSYS  finite-element  solutions  for  simulated  EM61  response  (with  expanded  time  scale)  to  steel 
objects.  Largest  response  is  for  vertical,  cylinder-like  object  (prolate  longitudinal)  because  the  greatest  amount  of 
steel  is  in  line  with  the  primary  field.  A  plate-like  object  on-end  (oblate  transverse)  similarly  has  a  large  response. 
Flat-lying  steel  plate  (prolate  transverse)  has  smallest  relative  response  because  of  unfavorable  orientation. 
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Figure  3.3.  Eddy  currents  in  81-mm  projectile  at  0.6  ms  (left)  and  9.6  ms  (right)  after  transmitter  turn-off.  Note 
strong  early  time  currents  in  tail-fin  shroud  and  late-time  concentration  in  tail-boom  threads.  The  latter  occurs 
because  currents  are  unable  to  diffuse  into  the  steel  tail-boom  screw  and  are  therefore  supported  to  longer  times  than 
expected  for  a  simple  ring. 


3.3  Dipole  Model 

Our  dipole  model  extends  earlier  work  (. McNeill  and  Bosnar,  1996;  Barrow  and  Nelson,  1999, 
2000;  Pasion  and  Oldenburg,  2001;  Machines  et  al.,  2002)  of  relatively  simple  models  for  UXO 
characterization.  The  target  is  treated  as  three  orthogonal  dipoles  with  one,  two,  or  three 
independent  axes  (a  sphere  has  one  independent  axis  and  an  axisymmetric  object  has  two).  The 
incident  magnetic  field  from  a  square  transmitter  loop  {Das  et  al.,  1990)  is  geometrically 
projected  upon  each  of  these  axes  depending  on  target  position  and  orientation.  The  time-domain 
EMF  V  in  each  axis  of  the  target  is  assumed  to  respond  as 

V  =  p{t+S)~rQat  (3.3) 

where  a,  /3,  y,  and  8  vary  for  each  independent  axis.  The  dipole  secondary  fields  from  each  target 
axis  are  evaluated  at  the  receiver  to  form  the  complete  response.  A  spherical  target  measured 
with  a  single-gate  system  (like  the  EM61)  requires  only  4  parameters  (/?i,  xyz  position)  whereas 
a  fully  time-dependent  triaxial  model  requires  18  parameters  («i23,  Pm,  Ym,  8m,  xyz  position, 
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azimuth,  inclination,  and  roll).  Grimm  and  Sprott  (2002)  applied  a  triaxial  model  in  ft  alone  (9 
parameters)  to  design  efficient  multicomponent,  multiple-sensor  platforms  for  UXO 
discrimination  (see  Wold  et  al.,  2002,  for  implementation).  Physically  the  ft  parameters  are 
proportional  to  the  time  derivative  of  the  magnetic  polarizability  along  each  principal  axis.  The 
a  parameters  are  the  “late-time”  exponential-decay  constants,  whereas  the  y  parameters  describe 
the  approximate  power-law  decay  that  would  be  represented  analytically  by  a  sum  of  higher- 
order  decay  constants  that  dominate  “intermediate  time.”  The  8  parameters,  introduced  by 
Pasion  and  Oldenburg  (2001),  are  effectively  empirical  correction  factors  to  improve  the  fit 
quality  (note  that  these  authors  actually  use  the  symbol  a  for  this  term).  In  contrast  to  the  other 
parameters,  clear  patterns  in  the  &  could  not  always  be  discerned  and  in  some  cases  the 
improvement  in  lit  caused  degradation  in  interpretation  of  ft,  a,  and/or  y.  The  (5-parameters  were 
therefore  treated  as  an  option  and  a  full  set  of  inversions  was  performed  both  with  and  without 
these  parameters. 

The  inverse  problem  for  target  shape  and  orientation  is  solved  through  a  generalized  inverse 
using  singular-value  decomposition  (e.g.,  Menke,  1984).  Direct  solutions  for  all  parameters  was 
generally  not  as  effective  as  the  following  multi-step  “recipe:” 

1.  Solve  for  the  best-fitting  sphere  (one  independent  axis  only)  in  all  relevant  parameters 
except  8. 

2.  Copy  these  estimates  to  the  other  axes  and  re-solve  with  3  independent  axes,  including  8  if 
desired. 

3.  Exchange  all  parameters  on  each  axis  so  that  the  largest  ft  value  is  the  1-axis  and  the 
smallest  ft  is  the  3-axis.  Re-solve  for  orientation,  holding  all  other  parameters  fixed. 

The  time-domain  response  decreases  sharply  with  time  and  also  varies  in  strength  with  spatial 
component  (Z  is  typically  larger  than  X  or  Y).  In  order  to  balance  all  of  the  data  in  the  least- 
squares  inversion,  time-  and  component-equalization  were  implemented  by  normalizing  each 
time  channel  and  each  component  by  its  respective  RMS  value  over  all  the  data  collected  for  that 
target.  Component  equalization  was  always  applied,  effectively  nulling  the  overall  Z-gain 
problem  described  above.  However,  this  approach  does  not  correct  any  intra-anomaly  variations. 
Time  equalization  was  treated  as  a  free  parameter;  a  full  set  of  inversions  was  performed  with 
and  without  time  equalization. 
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4.  ANALYSIS  OF  EM61-3D  TEST-SITE  DATA 


4.1  Instrumentation 

In  geological  EM  exploration,  the  optimum  approach  to  best  describing  target  size,  shape,  and 
orientation  is  to  use  multiple  time  gates  and  all  three  spatial  components  of  the  secondary  field 
(e.g.,  Kaufman  and  Keller,  1985).  The  Geonics  EM61-3D  is  a  multicomponent,  multichannel 
pulsed-EM  induction  sensor  designed  to  bring  these  capabilities  to  UXO  discrimination 
(Geonics,  Ltd.,  1996).  Previous  work  with  the  EM61-3D  and  related  prototypes  are  described  by 
McNeill  and  Bosnar  (1996)  and  Grimm  et  al.  (1997).  The  more-recent  Geonics  EM63  measures 
a  comparable  time  range  but  eliminates  recording  the  horizontal  spatial  components  in  favor  of 
higher  signal-to-noise  for  the  vertical  component.  The  EM61-II  also  only  measures  the  vertical 
component  in  a  more  restricted  number  of  time  gates.  Because  of  numerous  differences  in 
instrument  performance,  the  relative  merits  of  various  combinations  of  time  gates  and  spatial 
components  were  assessed  by  progressive  decimation  of  EM61-3D  data  rather  than  direct 
comparison  of  different  instruments. 

The  EM61-3D  consists  of  a  1-m  diameter  square  transmitter  coil  in  the  horizontal  plane  and  3 
orthogonal  receiver  coils  averaging  0.44-m  diameter  (Figure  4.1).  The  Geonics  PROTEM  is 
used  as  a  receiver.  The  transmitted  wavefonn  is  bipolar  quasi-rectangular  at  50%  duty.  The 
instrument  is  designed  to  operate  at  repetition  frequencies  of  7.5  and  30  Hz;  higher-quality  data 
was  obtained  at  the  lower  frequency  and  therefore  it  was  used  exclusively.  This  relatively  low 
pulse  rate,  and  perhaps  other  instrument  characteristics,  further  required  that  data  were  acquired 
statically,  i.e.,  at  a  series  of  fixed  stations,  rather  than  while  moving  (see  Grimm  et  al.,  1997). 
The  EMF  induced  in  each  receiver  coil  is  recorded  in  20  geometrically  spaced  gates  spanning 
0.35  to  27.9  ms  for  the  7.5-Hz  base  frequency. 

The  EM61-3D  shows  a  strong  time-dependent  response  even  when  there  is  no  metallic  target 
present.  The  Z-component  response  of  worst-case  conductive  ground  (10  Q-m)  is  <2  nV/m  at 
0.35  ms,  the  first  EM61-3D  time  gate  (TEMIX  XL).  However,  the  equivalent  Blossom-Point 
EM61-3D  data  (see  below)  show  target- free  responses  of  50-80  nV/nT.  Therefore  the  effect 
must  be  an  instrument  artifact,  which  was  approximately  removed  by  subtracting  the  mean  value 
around  the  periphery  of  each  measurement  grid  at  each  time  gate. 
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A  second  apparent  artifact  regarding  the  EM61-3D  concerns  the  relative  gain  of  the  Z- 
component  compared  to  X-  and  Y-components:  the  Z-component  appeared  to  be  too  small  in  the 
raw  data.  The  Z-component  of  the  analytical  response  of  a  sphere  ( Grant  and  West,  1965)  is  a 
factor  of  ~3  larger  than  the  EM61-3D  response  of  the  Blossom-Point  shotput,  when  compared  to 
their  respective  XY  components.  This  factor  was  therefore  applied  as  a  gain  to  the  EM61-3D  Z- 
component;  it  must  be  emphasized,  however,  that  this  value  is  an  average  for  measurements  of 
just  one  target;  it  may  not  be  constant  within  a  target  or  among  different  targets. 

4.2  Test  Objects  and  Measurements 

EM61-3D  measurements  were  made  at  the  test  range  established  by  the  Naval  Research 
Laboratory  (NRL)  at  Blossom  Point,  MD  (. Nelson  et  al.,  2000).  50-cal  and  20-mm  projectiles 
and  targets  deeper  than  0.5  m  were  excluded  for  potentially  low  signal-to-noise,  yielding  a  data 
set  consisting  of  25  unique  objects  in  45  different  states  of  depth  and  orientation  (i.e.,  typically 
two  instances  of  each  object  in  a  different  depth  and/or  orientation).  These  objects  included  60- 
and  81  -mm  projectiles,  a  Mk-23  practice  bomb,  a  16-lb  shotput,  steel  cylinders  and  rectangular 
plates  with  different  aspect  ratios,  and  numerous  clutter  items  (Table  4.1). 

Each  target  was  measured  at  73  discrete  stations  in  a  nonlinear  “starburst”  pattern  (Figure 
4.2),  that  is,  a  radial  pattern  in  12  directions  (30°  azimuth  separations)  at  distances  of  0.15,  0.3, 
0.5,  0.85,  1.5,  and  2.7  m,  in  addition  to  (0,0).  This  pattern  allowed  recovery  of  detailed  anomaly 
shapes  over  the  target  but  also  included  low-resolution  mapping  on  the  periphery  of  the  anomaly. 
Debiasing  and  gain  corrections  were  applied  as  discussed  above. 

4.3  Classification 

UXO  discrimination  or  classification  has  proved  to  be  more  effective  upon  model-derived 
parameters  rather  than  the  geometry  or  energy  of  the  geophysical  anomaly  itself  {Hart  et  al., 
2000;  Collins  et  al.,  2001).  Here  the  intrinsic  parameters  derived  from  the  dipole  model  (those 
that  are  properties  of  the  target  itself)  are  used  and  the  extrinsic  parameters  (derived  position  and 
orientation)  are  disregarded.  Direct  discrimination  upon  all  of  the  intrinsic  parameters  of  the 
time-dependent  triaxial  dipole  is  a  12-dimensional  problem:  although  classifiers  may  perfonn 
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well  on  a  training  data  such  as  that  used  here,  generalization  to  new  data  can  be  poor  because  the 
classifier  is  making  direct  comparisons  against  many  individual  parameters  for  each  training 
example.  Two  approaches  to  reducing  the  parameter  space  are  nondimensionalization  and 
selection  of  more  specific  discrimination  metrics  derived  from  the  individual  model  parameters. 

Nondimensionalization  is  implemented  by  dividing  each  parameter  by  a  reference  value. 
Recall  that  the  polarizability  derivatives  were  defined  above  so  A  >  A  >  /%;  adopting  the  1-axis 
as  the  reference,  the  three  J3s  are  expressed  as  two  parameters  /f, / A  <  A/A  ^  1.  Other 
parameters  are  expressed  similarly  with  respect  to  the  1-axis,  although  there  are  no  inequality 
constraints.  This  nondimensionalization  now  implies  that  classification  will  be  more  weighted 
by  shape — a  relative  measure — and  less  by  the  model  parameters  themselves,  which  indicate 
both  shape  and  absolute  size.  Although  there  may  be  less  difference  between,  say,  40-mm  and 
155-mm  projectiles  using  these  metrics,  the  goal  is  to  achieve  better  discrimination  of  long, 
axisymmetric  objects  relative  to  more  irregularly  shaped  scrap,  regardless  of  size. 

Indeed,  near-axisymmetry  is  a  defining  shape  characteristic  of  most  ordnance  and  is  adopted 
in  a  weak  and  strong  form  in  two  classification  algorithms.  The  weak  form  classifies  objects  by 
their  axisymmetry  but  then  simply  measures  how  well  the  shape  factors  above  are  grouped  by 
class.  The  strong  form  uses  the  same  classification  but  explicitly  measures  the  axisymmetry  of 
the  object  by  the  normalized  difference  in  shape  factors,  e.g., 

Ax,  =  l-IA/A-A/AIWA+A/A)  (4.  l) 

and  similarly  for  other  parameters  as  desired.  This  is  an  implementation  of  the  second  approach 
above  to  parameter-space  reduction.  The  minimum  value  of  the  axisymmetry  is  0  (when  either 
of  the  2-  or  3-axes  is  zero)  and  the  maximum  is  1  (when  the  values  for  the  2-  and  3-axes  are 
equal). 

The  Blossom-Point  test  objects  were  classified  as  “ordnance-like”  where  axisymmetric  about 
a  long  axis,  i.e.,  roughly  cylindrical  or  spherical  shapes.  The  10  distinct  ordnance-like  objects 
consisted  of  the  steel  cylinders,  projectiles,  practice  bomb,  and  shotput.  The  presence  of  tail  fins 
that  formally  break  true  axisymmetry  was  disregarded.  These  ordnance-like  objects  were 
measured  in  a  total  of  23  depth-orientation  states.  Other  items  that  were  triaxial  or  axisymmetric 
about  a  short  axis  were  classified  as  “scrap-like.”  The  15  scrap-like  objects  consisted  of  the 
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rectangular  steel  plates  and  clutter  items  described  above  and  were  measured  in  a  total  of  22 
depth-orientation  states. 

The  weak  form  of  the  axisymmetry  discriminator  was  implemented  by  a  generalized 
regression  neural  network  (GRNN;  Wasserman,  1993).  This  is  a  radial  basis-function  method 
that  requires  no  iterative  training  and  approximates  any  arbitrary  function  between  input  and 
output  vectors.  The  input  vectors  consist  of  up  to  8  shape  factors  described  above  (2  for  each  of 
the  4  parameter  groups  in  Eqn.  3.3).  As  a  classifier,  the  output  vectors  are  taken  to  be  the  class 
probabilities  and  the  GRNN  approximates  a  probabilistic  neural  network  (PNN;  Wasserman, 
1993;  Hart  et  al.,  2000).  In  a  two-class  system,  the  PNN  output  vectors  provided  for  training  are 
usually  taken  as  0  or  1;  for  the  GRNN  nonzero  values  are  preferred  numerically,  so  the  target 
outputs  are  assigned  1  or  2.  The  simplest  GRNN  or  PNN  has  a  single  free  parameter  u  which  is 
the  standard  deviation  of  the  gaussian  basis  function  or  kernel.  When  a  is  small,  the  GRNN 
extrapolates  and  the  PNN  classifies  based  on  the  nearest  neighbor  in  the  set  of  input  vectors. 
When  cr  is  large,  the  PNN  approximates  a  Bayesian  classifier  (i.e.,  it  yields  the  mode  of  the 
training  values)  but  GRNN  produces  the  (nonintegral)  mean  of  the  training  values.  The  network 
performance  is  assessed  using  leave-one-out  cross  training.  In  this  algorithm,  the  network  is 
asked  to  predict  the  output  vector  at  one  withheld  data  point.  Now  the  test  datum  could  be  most 
simply  classified  according  to  its  rounded  output  vector  (e.g.,  >1.5  is  predicted  ordnance-like  for 
the  GRNN)  but  in  practice  a  series  of  thresholds  is  tested.  At  each  threshold,  the  fraction  of 
correctly  classified  ordnance  (true  positives)  is  designated  the  Probability  of  Detection  (PD)  and 
the  fraction  of  scrap  incorrectly  classified  as  ordnance  (false  positives)  is  designated  the 
Probability  of  False  Alarm  (PFA).  The  ensemble  of  PFA  vs.  PD  as  a  function  of  output 
threshold  is  the  Receiver  Operating  Characteristic  (ROC)  curve. 

The  strong  form  of  axisymmetry  discrimination  uses  a  straightforward  cutoff  in  the 
axisymmetry  metric  (Eqn.  4.1).  Therefore  the  strong  cutoff  discriminator  uses  “hard”  weights 
whereas  the  weak  GRNN  discriminator  uses  “soft”  weights.  The  strong-form  algorithm  allows 
different  parameters  to  use  different  cutoffs  such  that  the  cutoff  surface  is  ellipsoidal,  i.e., 

R=X1[(l-^)/(l-c1)]2  (4.2) 
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where  R  is  the  normalized  distance  from  100%  axisymmetry  for  a  test  datum  (R<  I  is  classified  as 
ordnance),  A\  is  the  ith  of  up  to  4  axisymmetry  measures  for  that  datum  (e.g.,  /?axi)  and  0  <  c\  <  1 
is  the  cutoff  value  for  the  ith  parameter.  The  cutoff  is  swept  independently  for  each  of  the 
parameters;  the  fraction  of  ordnance  and  the  fraction  of  scrap  for  each  multiparameter  cutoff 
value  at  R  <  1  are  accumulated  into  PD  and  PFA,  respectively,  and  a  sort  is  subsequently 
required  to  order  PD  and  PFA  for  a  ROC  curve. 

4.4  Results 

4.4. 1  Sample  Parameter  Inversion 

Sample  results  are  shown  in  Figures  4.3-4.6  and  Table  4.2.  The  target  is  a  3”  diameter  by 
12”  long  stainless  steel  cylinder  buried  at  0.5-m  below  ground  surface,  lying  horizontally  and 
pointing  to  magnetic  north.  The  EM61-3D  data  were  interpolated  to  a  regular  grid  for  display 
but  analysis  is  performed  on  the  raw  data.  Note  the  distinct  character  of  the  X,  Y,  and  Z  spatial 
components.  The  sample  inversion  (Table  4.2)  used  the  full  data  set  under  spatial  and  temporal 
equalization  and  included  the  8  parameters.  The  overall  goodness-of-fit  is  96%.  The  ft  and  y 
parameters  for  the  minor  axes  are  79%  and  99%  axisymmetric  (Eqn.  4.1),  respectively,  but  the  a 
and  ^parameters  for  the  minor  axes  are  0%  and  43%  axisymmetry,  respectively.  The  zero  value 
of  axisymmetry  in  a  follows  from  a  zero  value  in  a  minor  axis,  which  in  turn  was  an  explicit 
movement  and  not  a  lack  of  change  from  initialization.  The  ratio  of  the  major  [4  to  the  geometric 
mean  of  the  minor  values  is  3.4:1,  which  somewhat  underestimates  both  the  aspect  ratio  of  the 
cylinder  (4:1)  and  the  DC  polarizability  ratio  (5.8:1)  for  a  magnetically  permeable  ellipsoid  with 
similar  dimensions  to  the  cylinder  ( Das  et  al.,  1990).  Recall,  however,  that  the  J3s  are  time 
derivatives  of  the  polarizability  and  are  also  dimensionally  influenced  by  values  of  y.  The 
recovered  position  and  orientation  are  in  excellent  agreement  with  the  true  values;  the  largest 
discrepancy  is  depth,  which  is  too  large  by  0.14  m  or  16%  of  the  sensor-target  distance. 
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4.4.2  General  Solutions 


Several  inversions  were  performed  for  each  of  the  45  buried  targets,  switching  by  system, 
parameterization,  and  equalization.  In  the  first,  EM61-3D  data  were  decimated  to  emulate 
systems  acquiring  less  spatial  and/or  temporal  data.  The  systems  were  denoted  as  follows: 

A  -  Z-component  receiver,  single  time  gate  centered  at  0.6  ms  (EM61  simulant). 

B  -  XYZ  receiver,  single  time  gate. 

C  -  Z  receiver,  20  time  gates  (EM63  simulant). 

D  -  XYZ  receiver,  20  time  gates  (full  EM61-3D). 

E  -  XY  receiver,  1  time  gate. 

F  -  XY  receiver,  20  time  gates. 

The  second  and  third  switches  were  whether  or  not  the  8  parameter  and  temporal  equalization 
were  selected.  In  all,  these  variations  led  to  675  individual  inversions. 

Results  are  grouped  by  instrument  configuration  and  discrimination  approach  (Table  4.3). 
The  column  headings  are  configuration  given  above.  The  row  headings  describe  the  model 
parameters  and  specific  method  used  for  discrimination.  “Ratios”  and  “LogRatios”  indicate  the 
weak-form  (GRNN)  discrimination,  whereas  “Axisymm”  denotes  the  strong-form  (axisymmetry 
cutoff)  discrimination.  The  average  PFA  over  all  inversions  at  91%  PD  is  reported  as  the 
principal  measure  of  discrimination  perfonnance  (This  PD  corresponds  to  21/23  correctly 
classified  ordnance-like  objects).  Configurations  A,  B,  and  E  cannot  estimate  the  temporal 
parameters  a ,  y,  or  8  and  so  these  entries  are  blank.  Good  performance  (highlighted  as  PFA  < 
50%)  largely  lies  within  configurations  D  and  F,  i.e.,  the  systems  with  multiple  time  channels 
recording  at  least  the  X  and  Y  field  components. 

The  value  of  8  and  time  equalization  were  assessed  by  the  average  changes  in  PFA  at  PD  = 
91%  by  alternately  incorporating  and  neglecting  these  factors.  As  both  deal  with  temporal 
effects,  they  apply  only  to  the  multichannel  configurations  C,  D,  and  F.  The  parameter  8 results 
in  a  small  but  quantifiable  improvement  in  discrimination  (Table  4.4),  whereas  time  equalization 
yields  a  major  increase  in  discrimination  ability  (Table  4.5)  because  the  temporal  parameters  are 
better  estimated.  Note  that  these  conclusions  are  based  on  typical  changes;  individual  cases  can 
be  adversely  affected. 
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4.4.3  Discrimination  Using  the  Power-Law  Slope 


The  final  results  (Table  4.6)  incorporate  both  8  and  time  equalization.  The  median  goodness 
of  fit  is  93%,  with  no  difference  between  ordnance-  and  scrap-like  items.  Very  good 
discrimination  performance  (PFA  <  25%  at  91%  PD)  is  highlighted.  The  best  classification  (5% 
PFA  at  91%  PD)  was  achieved  using  the  full  EM61-3D  multicomponent  multichannel  data  and 
analyzing  the  ratios  of  y  alone  (Table  6,  “/-Ratios”).  The  optimum  value  of  the  GRNN 
nonnalized  standard  deviation  o  is  0.7,  indicating  that  the  neural-network  weighting  spans  a 
significant  portion  of  the  parameter  space  and  therefore  should  be  somewhat  generalizable.  The 
quality  of  the  classification  is  intuitively  obvious  from  the  tight  grouping  in  the  /  ratios  of 
ordnance-like  objects  (Figure  4.7).  This  grouping  implies  values  just  greater  than  unity  for  all 
three  /,  which  is  shown  in  expanded  fonn  in  Figure  4.8a.  The  geometric  mean  of  the  three 
values  </>  (Figure  4.8b)  emphasizes  the  relatively  narrow  range  (1-2)  of  this  parameter  for 
ordnance-like  objects.  There  is  much  more  variability  in  /  among  scrap-like  items  in  both  of 
these  plots.  By  simply  selecting  </>  =  1-2  in  Figure  4.8b,  PD  =  100%,  PFA  =  39%  is  realized. 

The  close  clustering  of  / for  ordnance-like  objects  is  due  to  their  ability  to  sustain  diagnostic 
eddy  currents  longer  than  scrap-like  objects.  McNeill  and  Bosnar  (1996)  describe  the  different 
stages  of  time-domain  electromagnetic  induction  in  a  compact  permeable  conductor.  “Early 
time,”  just  after  transmitter  turn-off,  is  characterized  by  a  near-constant  magnetic  field  in  the 
target  that  is  established  by  surface  eddy  currents  as  a  consequence  of  Faraday’s  Law.  As  eddy 

_l/9 

currents  diffuse  into  the  object’s  interior  the  magnetic  field  decays  as  t  ;  this  is  “intermediate 
time.”  Finally,  eddy  currents  penetrate  the  entire  body  during  “late  time,”  which  is  characterized 
by  an  exponential  decay  of  the  magnetic  field.  When  EMF  is  the  measured  quantity,  voltage 
decays  as  t  during  intermediate  time.  The  beginning  and  end  of  intermediate  time  are 
determined  by  the  size  and  shape  of  the  object,  but  the  power-law  slope  of  the  EMF  during 
intermediate  time  is  always  -3/2.  The  observed  range  of  /between  1  and  2  for  the  ordnance-like 
items  is  in  excellent  agreement  with  theory.  As  an  object  transitions  to  late  time,  the  recovered 
value  of  /  should  in  principle  approach  0,  but  in  practice  a  wide  range  can  be  inferred  due  to 
ambiguity  in  the  curve  fitting.  This  is  exactly  what  is  observed  for  the  scrap-like  objects. 
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Pasion  and  Oldenburg  (2001)  suggested  that  a  mean  ;/ (which  they  call  p)  in  excess  of  0.8  for 
axisymmetric  objects  is  diagnostic  of  magnetically  permeable  materials.  They  calculated  the 
effect  of  variable  permeability  upon  spheres  of  different  sizes,  but  the  smallest  sphere  they 
considered  had  a  volume  of  ~500  cm  .  This  is  approximately  the  median  volume  of  the 
ordnance-like  objects  analyzed  here,  but  the  scrap-like  objects  have  a  median  volume  of  ~100 
cm  .  Only  three  of  the  scrap-like  items  are  nonferrous.  Therefore  the  result  of  Pasion  and 
Oldenburg  (2001)  should  be  generalized  to  state  that  magnetically  penneable  objects  show  /  in 
the  specific  range  1-2  when  the  object  is  sufficiently  large  to  sustain  early-to-intermediate  time 
currents  through  the  recorded  time  range.  Nonpenneable  or  smaller  permeable  objects  show 
widely  varying  /indicating  that  the  power-law  slope  was  not  accurately  recovered. 

The  comparatively  large  size  discrepancy  between  the  ordnance-like  and  scrap-like  items 
analyzed  here  (Figure  4.8c)  implies  that  the  demonstrated  clear  discrimination  in  y  is  not  due  to 
shape  under  the  weak-fonn  axisymmetry  but  to  size.  While  this  in  itself  is  a  powerful  and  useful 
result,  a  discriminant  is  still  desired  where  the  ordnance  and  scrap  are  comparable  in  size.  The 
strong  form  of  axisymmetry  can  be  so  applied,  as  it  effectively  measures  the  distance  to  the 
diagonal  in  each  panel  of  Figure  4.7  without  regard  to  clustering.  Figure  4.9  shows  that  the 
axisymmetry  separation  for  ordnance-like  vs.  scrap-like  objects  is  very  good  in  y,  good  in  /?,  and 
poor  in  a  and  8.  The  classification  ability  in  /alone  (Table  4.6,  “/-Axisymm”)  is  PFA  =  32%  at 
PD  =  91%.  This  measure  should  provide  greater  generalization  to  new  data. 

4.4.4  Ranking  of  Configurations 

The  mean  values  within  each  column  in  Table  6  can  be  used  to  assign  the  following  overall 
ranking  to  the  different  configurations: 

1  -  XYZ  components,  multiple  time  channels  (D) 

2  -  XY  components,  multiple  time  channels  (F) 

3  -  XYZ  components,  single  time  channel  (B) 

4  -  Z  component,  multiple  time  channels  (C) 

5  -  Z  component,  single  time  channels  (A) 

6  -  XY  components,  single  time  channel  (E) 
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An  alternative  metric  that  used  the  four  best  overall  discriminators  for  the  multichannel 
configurations  {pay  log  ratios,  fly  log  ratios,  y log  ratios,  y  ratios)  yielded  the  same  ranking. 

The  full  multicomponent,  multichannel  data  set  clearly  provides  the  best  overall 
discrimination.  Multichannel  XY  was  the  second  best  and  much  better  than  multichannel  Z; 
indeed,  single-channel  XYZ  outperformed  the  latter.  However,  single-channel  XY  was  the 
worst-performing  configuration.  The  details  of  these  rankings  may  lie  in  the  overall  quality  of 
the  EM61-3D  data,  in  the  limited  number  and  variability  of  test  targets,  or  in  the  applicability  of 
the  dipole  model  or  granularity  in  the  inversion  procedures.  What  is  clear  is  that  there  is  a  large 
decrease  in  performance  from  the  XYZ  and  XY  multichannel  configurations  to  all  others, 
especially  when  the  four  best  discrimination  schemes  are  selected.  Therefore  both  multiple 
spatial  components  and  multiple  time  gates  must  be  measured  for  optimal  UXO  discrimination. 
Complete  ROC  curves  for  the  two  best  GRNN  discriminators  and  the  best  strong-form 
axisymmetry  classifier  are  compared  to  GRNN  results  for  a  single-component,  single-channel 
system  in  Figure  4.10.  Improvements  at  91%  PD  are  40-60%. 
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Figure  4.1.  Geonics  EM61-3D  at  field  test  site.  Horizontal,  square  transmitter  coil  contains  three  orthogonal, 
circular  receiver  coils.  Transmitter  and  receiver  electronics  are  to  left. 


Figure  4.2.  Nonlinear  grid  used  for  both  high  spatial  resolution  near  target  center  and  complete  coverage  distant 
from  target.  See  text  for  details. 
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Figure  4.3.  Interpolated  map  view  of  the  X-component  of  the  EM61-3D  response  in  time  gate  5  (0.8  ms)  for  a 
horizontal,  3”xl2”  steel  cylinder  buried  at  0.5-m  depth.  Grid  is  oriented  true  north  but  target  points  magnetic  north. 
Full  color  scale  (blue  to  red)  is  -5  to  +5  nV/m2. 


Figure  4.4.  As  Figure  4.3,Y-component,  same  color  scale  -5  to  +5  nV/m2. 
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Figure  4.5.  As  Figure  4.3,  Z-component.  Full  color  scale  is  0  to  20  nV/m2. 


Log  EM  Response 


Figure  4.6.  Temporal  response  of  target  in  previous  figures,  evaluated  at  x  =  y  =  -0.25  m.  Log  scale;  notch  in  Y- 
component  response  is  sign  reversal. 
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Figure  4.7.  Dimensionless  minor-axis  model  parameters  for  full  EM61-3D  data  (configuration  D)  including  both  8 
parameter  and  time  equalization.  Ordnance-like  objects  are  red  crosses  and  scrap-like  objects  are  blue  x’s.  Four 
scrap  objects  lie  outside  the  plotted  P  ranges.  Note  good  separation  of  ordnance  and  scrap  in  P  and  excellent 
separation  in  y. 
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3  Gammas 


Figure  4.8.  Discrimination  of  ordnance-like  objects  using  power-law  slope  (y)  of  TDEM  decay.  A:  breakout  of 
minimum,  intermediate,  and  maximum  ^-values  for  ordnance-like  (numbers  1-23,  red)  and  scrap-like  (numbers  24- 
45,  blue)  targets.  B:  Geometric  mean  of  3  y-valncs  <y>  for  ordnance-like  (magenta,  dark)  and  scrap-like  (cyan, 
light)  objects.  <y>  lies  in  range  1-2  for  ordnance-like  objects.  C :  <y>  vs.  object  volume  for  ordnance-like  (red 
crosses)  and  scrap-like  (blue  x’s)  objects.  Correlation  of  <y>  with  target  type  is  largely  due  to  size  in  this  data  set. 
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R  Axisymmetry  a  Axisymmetry 


Figure  4.9.  Histograms  of  minor-axis  model-parameter  axisymmetry  for  full  EM61-3D  data  (configuration  D) 
including  both  <5 parameter  and  time  equalization.  Ordnance-like  objects  are  magenta  (dark)  and  scrap-like  objects 
are  cyan  (light).  Note  very  good  axisymmetry  separation  in  y,  good  separation  in  f:L  and  poor  separation  in  a  and  8. 
Higher  axisymmetry  is  effectively  a  shorter  perpendicular  distance  to  the  dotted  diagonals  in  Figure  4.7.  Cutoffs  in 
multiparameter  axisymmetry  form  basis  of  “strong”  discrimination  method. 
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Figure  4.10.  Receiver  Operating  Characteristic  (ROC)  curves  for  selected  instrument  configurations  and 
discrimination  approaches.  “Weak”  form  discrimination  for  ^ratios  has  excellent  performance  but  is  likely  not 
sufficiently  general.  “Strong”  form  discrimination  in  y  axisymmetry  is  very  good  and  is  likely  to  lead  to  better 
generalization. 
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Table  4.1.  Modeled  Test-Range  Objects 


Nr. 

Description 

Depth, 

m 

Azimuth, 

deg 

Inclina¬ 
tion,  deg 

Nr. 

Description 

Depth, 

m 

Azimuth, 

deg 

Inclina¬ 
tion,  deg 

ORDNANCE-LIKE 

SCRAP-LIKE 

1 

8 1  -mm  proj . 

0.35 

0 

0 

24 

12”x3”x74”  plate 

0.15 

0 

0 

2 

6”xl.5”  cyl. 

0.1 

0 

90 

25 

12”x3”x74”  plate 

0.25 

0 

0 

3 

12”xl.5”  cyl. 

0.5 

0 

90 

26 

6”xl.5”x74”  plate 

0.05 

90 

0 

4 

8 1  -mm  proj . 

0.5 

0 

0 

27 

12”x3”x74”  plate 

0.15 

0 

90 

5 

12”xl.5”  cyl. 

0.35 

0 

0 

28 

12”x3”x74”  plate 

0.25 

90 

0 

6 

Mk-23  bomb 

0.25 

0 

0 

29 

12”x3”x74”  plate 

0.15 

90 

0 

7 

6”x3”  cyl. 

0.1 

0 

0 

30 

6”xl.5”x74”  plate 

0.08 

0 

90 

8 

6”x3”  cyl. 

0.1 

0 

90 

31 

box  fin 

0.35 

0 

0 

9 

Mk-23  bomb 

0.25 

0 

90 

32 

shovel  blade 

0.23 

0 

0 

10 

Mk-23  bomb 

0.15 

0 

0 

33 

flat  soda  can  (Al) 

0.03 

0 

0 

11 

8 1  -mm  proj . 

0.5 

0 

90 

34 

twisted  steel 

0.23 

0 

0 

12 

12”x3”  cyl. 

0.35 

0 

90 

35 

banding  wire 

0.1 

0 

0 

13 

12”x3”  cyl. 

0.35 

0 

0 

36 

barbed  wire 

0.1 

0 

0 

14 

12”x3”  cyl. 

0.5 

0 

0 

37 

4”x4”x74”  Al  plate 

0.05 

0 

0 

15 

12”x3”  cyl. 

0.5 

0 

90 

38 

box  fin 

0.1 

0 

0 

16 

6”xl.5”  cyl. 

0.2 

0 

0 

39 

rocket  fins 

0.1 

0 

0 

17 

6”xl.5”  cyl. 

0.2 

0 

90 

40 

8”  round  steel  plate 

0.15 

0 

0 

18 

8 1  -mm  proj . 

0.35 

0 

90 

41 

2”x8”x74”  Al  plate 

0.05 

0 

0 

19 

12”xl.5”  cyl. 

0.35 

0 

90 

42 

fence-post  cap 

0.08 

0 

0 

20 

3”xl.5”  cyl. 

0.05 

0 

0 

43 

8”x8”x74”  plate 

0.25 

90 

90 

21 

Mk-23  bomb 

0.15 

0 

90 

44 

4”x4”x74”  plate 

0.05 

90 

90 

22 

6”x3”  cyl. 

0.25 

0 

0 

45 

8”x8”x74”  plate 

0.25 

135 

90 

23 

16-lb  shotput 

0.25 

- 

- 

Azimuth  and  inclination  are  with  respect  to  longest  target  axis.  Azimuth  measured  east  from  north,  inclination  up 
from  horizontal. 


able  4.2.  Model  Solution  for  3”xl2”  Steel  Cylinder 


A  =  21,000 

/I  =  7600 

Pi  =  5000 

«i  =  0.042 

a2  =  0.0 

«3  =  0.005 

n  =  0.85 

y2  =  1 .46 

Yi  =  1 .42 

c>i  =  0.14 

<%  =  0.31 

^  =  0.085 

x  =  0.02 

y  =  0.01 

z  =  0.64 

az=  172 

inc  =  4 

roll  =  12 

Known  values  are  x  =  y  =  0;  z  =  0.5,  az  =  0  (or  180),  inc  =  0,  roll  =  arbitrary. 

Units  are  (3-nV-msVm2;  a-ms"1;  y-dimensionless;  5-ms;  x,  y,  z-m;  az,  inc,  roll-degrees;  angles  measured 
with  respect  to  /?raxis. 
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Table  4.3.  PFA  @  PD  =  91%  —  Average  Over  All  Results 


Discrimination 

A 

B 

Configuration 

C  D 

E 

F 

p-Axisymm 

82 

46 

73 

45 

96 

49 

p-Ratios 

68 

73 

69 

72 

91 

67 

p-LogRatios 

86 

64 

80 

64 

77 

56 

y-Axisymm 

— 

— 

83 

68 

— 

66 

y-Ratios 

— 

— 

66 

43 

— 

51 

y-Log  Ratios 

— 

— 

69 

44 

44 

Pa-Axisymm 

— 

— 

82 

56 

— 

67 

Pa-Ratios 

— 

— 

65 

69 

— 

66 

Pa-LogRatios 

— 

— 

80 

60 

— 

53 

Py- Axisymm 

— 

— 

81 

48 

— 

48 

Py-Ratios 

— 

— 

63 

49 

— 

47 

Py-LogRatios 

— 

— 

66 

41 

— 

32 

Pay-Axisymm 

— 

— 

82 

58 

— 

59 

p-Axisymm 

— 

— 

63 

44 

— 

47 

P-Ratios 

— 

— 

61 

47 

— 

42 

“Axisymm”  is  “strong”  form  discrimination;  “Ratios”  and  “LogRatios”  are  “weak”  forms. 


Good  performance  (<50%  PFA)  is  highlighted. 


Table  4.4.  Average  PFA  Differentials  Due  to  Delta  Parameter 


Solution 

A 

B 

Configuration 

C 

D 

E 

F 

p-Axisymm 

... 

... 

23 

5 

... 

-7 

p-Ratios 

— 

— 

-7 

-21 

— 

2 

P-LogRatios 

... 

... 

-18 

-18 

... 

-7 

y-Axisymm 

... 

... 

-11 

-18 

... 

-23 

y-Ratios 

— 

— 

-5 

9 

— 

2 

y-LogRatios 

... 

... 

-16 

7 

... 

2 

Pa-Axisymm 

... 

... 

18 

-16 

... 

-25 

Pa-Ratios 

— 

... 

7 

-16 

... 

-5 

Pa-LogRatios 

... 

... 

18 

-16 

... 

-2 

Py-Axisymm 

... 

... 

-7 

-5 

... 

-14 

Py-Ratios 

... 

... 

-2 

-16 

... 

-7 

Py-LogRatios 

... 

... 

0 

-14 

... 

-14 

Pay-Axisymm 

... 

... 

18 

-11 

... 

-18 

Pay-Ratios 

... 

... 

7 

-2 

... 

-2 

Pay-LogRatios 

... 

... 

-14 

-16 

... 

-16 

Median  changes  to  configurations  C,  D,  and  E  are  -2,  -16,  and  -7,  respectively. 

Other  configurations  cannot  model  temporal  parameters  and  therefore  are  not  applicable. 
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Table  4.5.  Average  PFA  Differentials  Due  to  Time  Equalization 


Solution 

A 

B 

Configuration 

C 

D 

E 

F 

p-Axisymm 

... 

... 

-27 

0 

... 

20 

p-Ratios 

— 

— 

-20 

-7 

— 

-7 

p-LogRatios 

... 

... 

-11 

-5 

... 

-25 

y-Axisymm 

... 

... 

-16 

-46 

... 

-14 

y-Ratios 

— 

— 

-39 

-50 

— 

-48 

y-LogRatios 

... 

... 

-36 

-43 

... 

-43 

Pa-Axisymm 

... 

... 

0 

18 

... 

23 

Pa-Ratios 

— 

— 

-32 

-16 

— 

0 

Pa-LogRatios 

... 

... 

-25 

-7 

... 

-30 

Py-Axisymm 

... 

... 

-18 

0 

... 

5 

Py-Ratios 

... 

... 

-25 

-30 

... 

-25 

Py-LogRatios 

... 

... 

-43 

-27 

... 

-32 

Pay-Axisymm 

... 

... 

0 

9 

... 

41 

Pay-Ratios 

... 

... 

-34 

-43 

... 

-25 

Pay-LogRatios 

— 

— 

-48 

-30 

— 

-43 

Median  changes  to  configurations  C,  D,  and  E  are  -25,  -16,  and  -25,  respectively. 

Other  configurations  cannot  model  temporal  parameters  and  therefore  are  not  applicable. 


Table  4.6.  PFA  @  PD  =  91%.  Final  Results  Using  Delta  Parameter  and  Time  Equalization. 


Solution 

A 

B 

Configuration 

C  D 

E 

F 

p-Axisymm 

82 

46 

86 

55 

96 

55 

p-Ratios 

68 

73 

77 

59 

91 

77 

P-LogRatios 

86 

64 

77 

50 

77 

46 

y-Axisymm 

— 

— 

100 

32 

— 

55 

y-Ratios 

— 

— 

64 

5 

— 

27 

y-LogRatios 

— 

— 

68 

9 

23 

Pa-Axisymm 

... 

... 

100 

100 

... 

100 

Pa-Ratios 

... 

... 

64 

55 

... 

68 

Pa-LogRatios 

... 

... 

77 

59 

... 

46 

Py-Axisymm 

... 

... 

68 

50 

... 

41 

Py-Ratios 

— 

— 

73 

23 

— 

46 

Py-LogRatios 

— 

— 

77 

9 

23 

Pay-Axisymm 

— 

— 

100 

100 

— 

100 

Pay-Ratios 

— 

— 

64 

23 

— 

46 

Pay-LogRatios 

— 

— 

46 

23 

18 

High  performance  (<25%  PFA)  is  highlighted.  100%  PFA  indicates  ROC  curve  did  not  achieve  91%  PD. 
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5.  CONCLUSIONS 


Of  the  three  physics-based  models  investigated,  only  the  simplest — the  dipole — is  presently 
suitable  for  UXO  studies.  We  extended  this  model  from  the  triaxial,  time-independent  model  of 
Barrow  and  Nelson  (2001)  and  the  axisymmetric,  time-dependent  model  of  Pasion  and 
Oldenburg  (2001)  to  a  full  triaxial,  time-dependent  representation.  This  model  was  used  to  infer 
target  properties  using  multicomponent,  multichannel  time-domain  electromagnetic  signatures  of 
45  seeded  objects  were  obtained  using  the  Geonics  EM61-3D.  Instrument  limitations  required 
static  data  acquisition;  a  nonlinear  grid  was  used  to  assure  both  high  spatial  density  directly  over 
the  targets  and  adequately  broad  coverage.  Neither  of  these  constraints  is  applicable  to  practical 
UXO  field  work.  However,  they  demonstrate  optimum  discrimination  ability  and  use  a 
prototype  instrument  that  provides  an  internally  consistent  basis  for  assessing  the  value  of 
different  combinations  of  spatial  components  and  multiple  time  gates.  The  latter  was 
implemented  by  decimating  the  data  by  component  and  channel  to  provide  several  alternative 
configurations  acquiring  one,  two,  or  three  spatial  components  and  one  or  multiple  time  gates. 

The  seeded  targets  were  divided  into  “ordnance-like”  and  “scrap-like;”  the  former  were 
axisymmetric  about  a  long  axis  whereas  the  latter  were  not.  The  triaxial  time-dependent  dipole 
model  was  used  to  estimate  up  to  18  parameters  regarding  target  size,  shape,  position,  and 
orientation.  When  nondimensionalized  by  the  values  along  the  object’s  apparent  longest  axis, 
there  are  8  parameters  relevant  to  target  size  and  shape.  A  generalized  regression  neural  network 
was  able  to  separate  the  ordnance-like  items  with  a  91%  true  positive  and  5%  false-positive  rate 
using  the  full  multicomponent,  multichannel  data  set.  However,  this  was  largely  due  to  the 
larger  sizes  of  the  ordnance-like  items  compared  to  the  scrap-like  objects:  the  voltage  decay 
approximately  follows  the  -3/2  log-log  slope  that  is  theoretically  expected  for  larger  targets  that 
are  able  to  extend  “intermediate  time”  through  the  measured  time  range.  Further  work  on  objects 
of  comparable  volume  but  different  shape  is  necessary  to  determine  whether  the  scrap-like 
objects  still  move  through  intermediate  time  more  quickly  because  they  generally  have  one  very 
thin  dimension.  An  alternative  metric  for  the  power-law  slope  that  depends  solely  on 
axisymmetry — without  regard  to  size — still  yielded  good  discrimination:  32%  false-positives  at 
91%  true  positives  in  the  ordnance-like  objects. 
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Multiple  spatial  components  and  multiple  time  channels  were  typically  able  reduce  the  false¬ 
positive  rate  by  40%  or  more  over  multiple-component,  single-channel  or  single-component, 
multiple-channel  decimated  data  sets,  demonstrating  the  improved  value  of  such  systems  for 
UXO  discrimination. 
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